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Abstract 

Two new test statistics are introduced to test the null hypotheses that the sampling distri- 
bution has an increasing hazard rate on a specified interval [0,a]. These statistics are empirical 
Li-type distances between the isotonic estimates, which use the monotonicity constraint, and 
either the empirical distribution function or the empirical cumulative hazard. They measure the 
excursions of the empirical estimates with respect to the isotonic estimates, due to local non- 
monotonicity. Asymptotic normality of the test statistics, if the hazard is strictly increasing on 
[0, a], is established under mild conditions. This is done by first approximating the global empir- 
ical distance by an distance with respect to the underlying distribution function. The resulting 
integral is treated as sum of increasingly many local integrals to which a CLT can be applied. 
The behavior of the local integrals is determined by a canonical process: the difference between 
the stochastic process x ^ W{x) + where W is standard two-sided Brownian Motion, and 
its greatest convex minorant. 

1 Introduction 

One way of characterizing a distribution of an absolutely continuous random variable X that is 
particularly useful in reliability theory and survival analysis, is by its hazard rate h^. Suppose X 
models the failure time of a certain device. The interpretation of the hazard rate is that for small 
e > 0, eho{x) reflects the probability of failure of the device in the time interval {x,x + e] given 
the device was still unimpaired at time x (assuming /iq is continuous at x). Put differently, hQ(x) 
represents the level of instantaneous risk of failure of the device at time x, given it still works at 
time X. A high value reflects high risk, a low value low risk. Lifetimes of devices that are subject 
to aging can be described by distributions with increasing hazard rate. Locally decreasing hazard 
rates can be used to model life times of devices that become more reliable with age during a certain 
period of time. 

It is especially this clear interpretation of these qualitative properties of a hazard rate that makes 
this function a natural characteristic of a survival distribution. The problem of estimating a hazard 
rate nonparametrically under qualitative (or shape) restrictions gained attention in the sixties of 
the previous century (see Groeneboom and Jongbloed (2011b) and the references therein). Also 
the problem of testing the null hypothesis of constant hazard (exponentiality) against monotonicity 
of the hazard was studied intensively, see e.g. Proschan and Pyke (1967). Only quite recently 
another testing problem, with a "shape constraint" rather than parametric null hypothesis was 
studied. See also the discussion in companion paper Groeneboom and Jongbloed (2011a). 
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In this paper, wc consider the asymptotic distribution theory for two integral-type test statistics 
for the hypothesis that a hazard rate ho is monotone on an interval [0,a], for some known a > 0. 
We restrict ourselves to the increasing case; the case of locally decreasing hazard can be considered 
analogously. 

Based on an i.i.d. sample Xi, . . . , Xn from the distribution associated with Hq, the most natural 
nonparametric estimator for Hq without assuming anything on Hq, is the empirical cumulative 
hazard function given by 



^ _ / - log {1 - F„(x)} , X G [0, , 
" I oo, x> 



where F„ denotes the empirical distribution function based on Xi,X2, . ■ ■ , Xn. Under the assump- 
tion that Hq is convex on [0, a], the cumulative hazard can be estimated by the greatest convex 
minorant Hn of the empirical cumulative hazard function 1HI„ on the interval [0, a]. Using these two 
estimators, the following test statistic emerges: 

Tn= I {Unix-) - Hn{x)] dFn{x). (1.1) 
J[0,a] 

Note that this is the empirical Li-distance between the two mentioned estimators for the cumulative 
hazard function w.r.t. the empirical measure dFn, and that r„ > since Hn is a minorant of 1HI„. 

If Hq is concave on [0, a], both estimators for Hq will be close to i^o a-^d T„ will tend to be small 
(converge to zero a.s. for n oo). On the contrary, if Hq has a region in [0, a] where it is not 
increasing, M.^ will capture this "non-convexity" of Hq and converge to Hq on this region whereas 
Hn will converge to the convex minorant of Hq on [0, a]. Note that r„ = if and only if Hn coincides 
with the linear interpolation of the points (x(j), ]HI„(a:(j) — )) on the range of the data falling in [0, a]. 
One could say that r„ = if ]HI„ is 'as convex as it can be on [0,a]', being an increasing right 
continuous step function. This is the reason for taking BI„(x— ) instead of ]HI„(x) in (1.1). A similar 
reasoning can be held for another test statistic, 

Un= [ {F„(x-) - Fnix)} dFnix), where F„(x) = 1 - exp(-i?„(x)). (1.2) 

An advantage of this definition is that Un is less sensitive to possible problems that can occur with 
large values of 1HI„. 

The main result of this paper concerns the asymptotic distribution of T„ and Un- under certain 
assumptions 

„5/6 I _ ^^^1 ^ ^ ^2^J ^5/6 _ ^jj^y ^ ^ (q^ ^2 J ^ 3) 

where Tn is a modified version of r„, see Theorems 4.1 and 4.2. Here ajj^ and ajp^ are constants 
depending on /q. Results of a similar flavor were established in, e.g., Kulikov and Lopuhaa (2008) 
for the difference between the empirical distribution function and its concave majorant. 

The basic idea of the proof is to approximate the integral in the test statistic as sum of in- 
creasingly many local integrals, using the crucial localization Lemma 3.4, and to apply a Central 
Limit Theorem to the components that arise in this way. The behavior of the local integrals is 
determined by a canonical process, the difference between a Brownian motion with parabolic drift 
and its convex minorant. Relevant properties of this process are derived in section 2. In section 
3, a statistic related to r„ (where the integral is taken w.r.t. Fq rather than F„) is closely approx- 
imated by an integral involving the independent increments of Brownian motion. Moreover, the 
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resulting integral is represented as a sum of local integrals using a "big blocks separated by small 
blocks" construction as introduced in Rosenblatt (1956). The local integrals over the big blocks 
reduce to the processes considered in section 2. Finally, because the local integrals are based on 
the independent increments of a Brownian motion process, a CLT can be applied to obtain the first 
result in the spirit of (1.3). In section 4, the main results of the paper are established by showing 
that the differences between the integrals w.r.t. (iF„ and cIFq are sufficiently small. 

2 Asymptotic local problem 

Consider the process 

X ^ V{x) = W{x) + x'^,xeR (2.4) 
with W standard two-sided Brownian motion on M. Then, for c > 0, define the functional Qc by 

Qc= [ {V{x) - C{x)} dx, (2.5) 
Jo 

where C is the greatest convex minorant of F on M. For a picture of the process V and its greatest 
convex minorant, restricted to the interval [—2, 2], see Figure 1. We have the following result. 

Theorem 2.1 

c~^/^ {Qc - cE\C{0)\} A A^(0, (7^), c ^ oo, 
where C (0) is the value of the greatest convex minorant C of the process V at zero, and 

a^ = 2 covar(-(7(0), V{x) - C{x)) dx. 
Jo 

All moments of c~^/^ {Qc — cE\C{0)\} exist and (in particular) the fourth moment is uniformly 
bounded in c and converges to the fourth moment of the normal N{0, cr^) distribution, as c oo. 

In the proof we will use the following lemma, which is proved in the appendix. 

Lemma 2.1 For the process V defined in (2.4), there exist positive constants c and c' such that 
for all u >0 



P min V{x) < < ce-" " . 

\x^[-u,u\ J 

Proof of Theorem 2.1. It follows from the results in Groeneboom (1989) that the process 

V{x) - C{x), xeR, (2.6) 

is stationary. In fact, the process touches zero at changes of slope of C and behaves between these 
touches of zero as an excursion of a Brownian motion path above a parabola of the form 

(j){x) = s — {x — a)2, X G M, 

where is a parabola touching two local minima of Brownian motion, and where the (random) 
values a and s depend on the Brownian motion path. Defining 

rk+l 



Dk= {V{x) - C{x)} dx, keZ, 

J h 
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Figure 1: The greatest convex minorant of W{x) + x^, restricted to [—2, 2]. 
we get a stationary sequence of random variables, and the stationarity of the process (2.6) yields: 

rk+l 



EDk= E{V{x) -C{x)} dx = E\C{0)\. 

Jk 

Moreover, all moments of exist. This follows from the fact that 

max {V{x) — C{x)} 
xe[o,i] 

has a distribution with tails which die out faster than exponentially. To see this, note that, \/u > 0, 

pi max|y(x) -C(x)} > mI < pi max V(x) > I + piminC(x) < -^mI 
Ue[o,i] J U-e[o,i] ^ J [x6iR ^ J 

< pi max W(x) > iM - 1 i + pi minF(x) < -^M \ 
{xe[o,i] ^ ' ^ J [xeR ^ J 

<\-[ e'S'^'dx + pi min W{x) < -^mX +F \ min V{x)<o\. (2.7) 
V vr J^M-i [xe[-u,u] J [x^[-u,u] J 

The first term on the right hand side is bounded by cexp{— c'M^/4} for some c, c' > 0. By 
Lemma 2.1 we have for the third term: 



P min V{x) < < ce 



X9:\~U,U\ 
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for constants c, d > 0. For the second term in (2.7), we get by Brownian scaling, 

min W(x) < -M/2 i = P i min W(ux) < -M/2 \ 
e[-u,u] J ta;e[-l,l] J 



= pj min u-^/'^W(ux) < -u-^/^M/2\ =¥1 max W(x) > u'^/'^M/l} 
Ue[-i,i] 'J \xe[-i,i] ' } 

< 2./^ r e-h' dx < 2./^^ exp {-iM^u} . 



Hence, taking u = v M in the second and third term in (2.7) and observing that the first term is 
of lower order, we obtain 

P{ max {V(x) - C(x)} > M} < cie-"^^^'^ (2.8) 

a;e[0,l] 

for constants ci, C2 > 0. 

Now let r(a) be defined by: 

T(a) = argmin^gK {W{x) + {x - af] . 

The (stationary) process a i— T(a) — a is studied in Groeneboom (1989) and it follows from the 
results, given there, that there exist positive constants ci and C2 such that 

|P(A n B) - P(A)P(B)| < 016-^=2™', 

for events A and B such that 

A^a {T{a) : a < 0} , 5 G a {r(a) : a > m} . 
This implies that there also exist positive constants c\ and C2 such that 

|P(A n S) - P(A)P(S) I < cie-^2m3 ^ ^2.9) 

for events ^ and B such that 

^ € o- {V{x) - C{x) : X < 0} , 5 G a {V{x) - C{x) : x > m} . 
So we can apply Theorem 18.5.3 in Ibragimow and Linnik (1971), p. 347, yielding that 

c-^/'{Qc-c£;|C(0)|} AiV(0,a2), 

where 

oo 

cr^ = var(L'o) + 2 ^ covar(Z)o, i^fe)- 
fe=i 

Using the stationarity of the process (2.6) again, we obtain 

fOO 



/•oo 

a^ = 2 covar(-C(0), V{x) - C{x)) dx. 
Jo 



'0 

The last statement of the theorem follows from (2.8) and (2.9). □ 
We will also need the following extension of Theorem 2.1. 
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Theorem 2.2 Let Cc be the greatest convex minorant on [0, c] of the process 

V{x), X e [0,c]. 

Note that Cc is not the restriction of C to [0, c], since C is globally defined on M, and Cc is the 
greatest convex minorant of the process V on [0,c], and only defined on [0, c]. 

(i) Let, for c> 4, the interval Lc be defined by 

Ic = [-\/c, c- ^/c\ . 

Then: 

c-^/^l^ {V{x)-Cc{x)} dx-E {F(x)-Cc(a;)}da;| AiV(0,a2),c^oo, (2.10) 

where is defined as in Theorem 2.1. 
(a) Relation (2.10) also holds if the interval I c is given by: 

Lc = [O, c - -s/c] or Lc = [s/c, c] . 

(Hi) For any choice of Lc in (i) or (ii), the fourth moment of 

c-^l'^ {/ ^^"^^ ~ '^"'^^^^ Ii ^^^^^ ~ ^"''^'^^ 

is uniformly bounded in c, and converges to the fourth moment of a normal N{0, a^) distri- 
bution, as c oo. 

Proof, (i). The probability that Cc is different from C on the interval Lc is less than or equal to 

fciexp 1-^20^/^} , 

for constants ki,k2 > 0. The proof of this is analogous to the proof of Lemma 3.4 in the next 
section. Hence, if Kc denotes the event that Cc ^ C on Ic, we get: 

E \Cc{x) - C{x)\ dx < 1^ E{V{x) - C{x)f dxj ^ F{Kc)^/^ = O (c^/^e-'"''^'^ , c ^ oo, 

for some A; > 0. Hence: 

1^ {^(2;) _ (^^(^^1 dx-E {V{x) - Cc{x)] dxj 

and the statement now follows. 

(ii). We can repeat the argument on the interval [0, ^/c\, and apply the argument used in (i) on the 
subinterval I'c = [c^^^, ^/c— c^^^] (but leaving Cc as it was defined in (i)). This yields: 
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implying: 

.-1/2 



1^ - Cc{x)} dx-E {V{x) - Cc{x)} dxj ^ 0, c ^ oo, 



Moreover, 

^1/2 / E \V{x) - Cc{x)\ dx = fc-V4) , c ^ oo. 
V[0,cV4] V / 



The statement now follows for the first choice of the interval Ic in (ii). For the second choice of Ic 

the argument is similar. 

(iii). Let Ic be as in (i). Then: 

c-'^E 1^ {V{x) - Cc{x)} dx-E {V{x) - Cc(x)} darj 

= c-'^E {V{x) - C{x)} dx-Ej^ {V{x) - C{x)} dxj + O (e"*^^'^') , 
for some k > 0, and the statement now follows from Theorem 2.1, (2.8) and (2.9) and the fact that 

f c y ^ 

If, for example, Ic = [0, c — y^, we write 

/ {V{x) - Cc{x)} dx-E f {V{x) - Ccix)} dx = Ac + 

Jlc Jio 

where 

dx 



Ac= [ {V{x) - Cc{x)} dx-E [ {V{x) - Cc{x)} 
Bc= I {V{x) - Cc{x)} dx-E I {V{x) - Cc{x)} dx. 

J k/c,c— \/cl J \\/c,c—^/c\ 



'[yfc,C-yfc\ J[yfc,C-^ 

Hence we get: 



Elf {V{x) - Cc{x)} dx-E I {V{x) - Cc{x)} dx 

[Jlc Jlc 

= c-'^EB^ + {4.EB^Ac + QEB^A^ + AEBcAl + EA^} . 

We have: 

c-'EAt=(^^Jc-'EAt = 0{c-'), 

and similarly, using the Cauchy-Schwarz inequality, 

c-^EBcAl = Ec-^/^Alc-'/^Bc < ,/Ec^,^/Ec^ = O (c'^/^) . 

Continuing in this way, we find that the only non- vanishing term is c~^EB^. The statement now 
follows from what we proved for Ic = [^/c, c — ^/c] . □ 

We finally also need the following extension of Theorem 2.2. 



Theorem 2.3 Let Fc, Gc and He be twice differentiable increasing functions on [0,c], with contin- 
uous derivatives fc, gc and he, respectively, satisfying 

Fc{x) = /c(0)x(l + o(l)), Gc{x) = 5c(0)x(l + o(l)), Hc{x) = ^h'M^^i^ + o{l)), c ^ oo, 

where the o(l) term is uniform in x. We assume that /c(0), gc{0), hc{0) and h'^{0) are positive and 
stay away from zero and oo, as c—^ oo, where h'^{0) denotes the right derivative of he at zero. Let 
Cc be the greatest convex minorant on [0, c] of the process 

Vcix) = Hc{x) + W{Gc{x)), X G [0,c]. 

Moreover, let Sc be defined by 

Sc{x) = Vc{x) - Cc{,x), X e%c]. 

Then: 

(i) Let, for c > 4, the interval Lc be defined by 

Ic = [Vc, C - Vc] • 

Then: 

c-'E [ Sc{x)dFc{x) ~ g^(")'"/^^(Q) ^|c(0)|, var fc'^/^ f s,{x) dFc{x)) ^ alc^oo, 
•lie {Ih'AG])'^ V Jic J 

(2.11) 



where 



2 5c(0)VVc(0)2 2 



and C and are defined as in Theorem 2.1. Moreover, the fourth moment of 

f {S^(^x) - ESc{x)} dFc{x) 

JIc 

is uniformly bounded, and satisfies: 

E (^-^/"^ {Sc{x) - ESc{x)} dFcix)^ ~ M^^^), c ^ oo, (2.13) 

where Mc^^ denotes the fourth moment of a normal N{0, a^) distribution. 

(ii) Relations (2.11) and (2.13) also hold if the interval Lc is given by: 

Lc = [0,c - y/c] , Lc= [s/c, c] or Lc = [0, c]. 

Proof. Since the proof proceeds along the hnes of the proofs of Theorems 2.1 and 2.2, we only pay 
attention to the new type of scaling which is present in the process 

X i-> |/i.c(0)ar^ + W{gciO)x), x G [0, c], 

which replaces the process 

X I— 7- V{x) = x^ + W{x), X G [0,c]. 
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Let a,b > 0. By Brownian scaling, the process 

X ax^ + W{bx), X e [0, c], (2.14) 
has the same distribution as the process 

I (|„2/3^-l/3^^\ ^(^2/3^,-1/3^) I ^ ^ ^ [q^ ^] (2.15) 

Hence, if Ca.b is the greatest convex minorant of the process given in (2.14) and Ca,b of the process 
given in (2.15) we get: 



r {ax^ + W{bx) - Ca,b{x)} /c(0) dx 
Jo 



'0 

V 



bfciO) 



' ' + W{u) - a'/^-^/^Ca,b (a-2/351/3^) I du 



a Jo 

where Cc is the greatest convex minorant of the process 

u^u^ + W{u), u e [0, a^/^b-^/^c 
Thus, for c — )• 00 

c~^e£ {ax^ + W{bx) - Ca,b{x)} /e(0) dx ~ ^!^^M^£;|C(0)|. 

Using that a = ^h'^iO), b = gdO), (2.11) foUows. Moreover, 
.-1/2 



var I c 



r {ax^ + iy(6x) - Ca,bix)} /e(0) 
JO y 

var / [u^ + W{u)-Cc{u)] du 



a^c Wo 



,4/3 



a' 



4/3 



-cr c — > 00. 



Taking a = ^^(0), b = gdO) now yields (2.12). □ 

3 Embedding and first central limit result 

In this section, a central limit result is established for the quantity 

£ {h„(x) - Hn{x) - fjLn} dFo{x) = u^/^ ^ {h„(x-) - ^„(x) - dFo{x), (3.16) 
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n 



where /u„ denotes a centering sequence to be specified below. This result is the first step to be taken 
in order to obtain the limit result for Tn defined in (1.1) (where the integral is taken with respect 
to dFn rather than cIFq). In order to derive the limiting distribution of (3.16), we first replace the 
process M.n{x) — Hn{x) by 

X ^ Ho{x) + ^^/'^""l , . - Hn{x), X e [0,a]. 



where En is the empirical process -^/nlF^ — Fq} and Hn is the greatest convex minorant of the 
process 

x^Ho{x) + n-^/^En{x)/{l-Fo{x)},xe[0,a]. (3.17) 

Next we use the strong approximation of the empirical process by a Brownian bridge yielding 
the approximation 

, , n-^/^BjFo(x)) ,^ , 
X ^ Ho{x) + — ;y X e [0,o], 

the process (3.17). This process is distributed as 

X ^ Ho{x) + n-'/^W (3^^^) ' ^ e [0> «]> (3-18) 

where W is standard Brownian motion on [0, 00). Next, the interval [0, a] is split up in so-called 
big blocks separated by small blocks. The local contributions to the integral over the big blocks 
can be treated using the results of section 2. 

The first lemma to be proved states a contraction property for convex minorants that will be 
used repeatedly in the sequel. It is related to Marshall's Lemma in the theory of isotonic regression. 

Lemma 3.1 Let f and g he hounded functions on an interval 7 C M and let Cf and Cg he their 
greatest convex minorants, respectively. Then: 

sup \Cf{x) - Cg{x)\ < sup \f{x) - g{x)\. 
xei xei 

Proof. Using that f > g — sup^^^j \f{u) — g{u)\ and that 5 > Cg by definition, it follows that 

f > Cg — supygj |/(^) ~ 5(^)1- Since the right hand side is convex, this means that it is a convex 
minorant of / on /. Hence, it lies below the greatest convex minorant C/ of / on I: 

Cf{x) > Cg{x) - sup |/(n) - g{u)l xEl. 

uei 

Since this inequality also holds with / and g interchanged, the result follows. □ 
We now consider the functional 

[ iUnix)- Hnix)} dFo{x) 
J [0,0] ^ 



X 



. Ho{x) - log (1 - ^.f"^!!, , J - Hn{x) \ dFoix), (3.19) 
[0,a] [ \ ^Al{l-Fo{x)} J J 

where = ^/n{¥n — Fq} is the empirical process. The following lemma enables us to dispense 
with the logarithms. 
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Lemma 3.2 Let Hn he the greatest convex minorant of the process 



X ^ Ho{x) + ^rf"^!!, ,v X e [0,a], 



where Fo{a) < 1. Then: 
(i) 



i 

J [0,a 



ln{x) - Hoix) - 



En{x) 



V^{l-Fo(x)} 



dFo(x) = Op (n-i) . 



(ii) 



I Hn{x)-Hn{x) dFo{x) = Op {n-^) 

J[0,a] 



Proof, (i). Let An denote the event 



EJx) 



sup , . 
ce[o,o] V'^ll - Eo{x)\ 

Then, by a well-known result in large deviation theory ("Chernoff's theorem"), we have 

P«) = 0(e— ), 

for a constant c > 0. If An occurs, we can expand the logarithm, which yields: 
and (i) now follows. 

(ii). This follows from Lemma 3.1 and the argument of the proof of (i). 



< 



□ 



n 



We shall prove below that 

V6 f I Ho{x) + ^,f"^t^ ^^ - Hn{x) - e\ Ho{x) + ^.^^"^^^ - Hn{x) } } dFo{x) 



'[o,a][ V^{l-Foix)} ^ I ' v^{l-Fo(x)} 

converges in distribution to a normal distribution, which, together with Lemma 3.2, implies that 



n^/e f Unix) - Hn{x) - e{ho{x) + ^,f"^!! 
Jlo,a] [ [ ^/n{l-Fo 



- Hnix) } } dFoix) 



(x)} 

converges to the same normal distribution. 
Remark 3.1 We avoid taking the expectation of 

Mn{x)-Hn{x), 

since ]HI„ is infinite with a positive (but vanishing) probability on [0,a], as is Hn- This happens 
when the empirical distribution function reaches the value 1 on [0, a]. 
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By Theorem 3 of Komlos, Major and Tusnady (1975) we can construct Brownian bridges 
on the same sample space as F„ such that 



Yn = sup 

xe[0,a] 



n^/^ \En{x) - Bn{Fo{x))\ 



2 V log n 

is a random variable with with EYn < C < oo for all n. Hence, for n>2: 



0<E sup n"^/2 

x&[0,a] 



En{x) Bn{Fo{x)) 



n(l - Fo(a)) 



l-Fo(a;) l-Fo(x) 
We now have the following result. 
Lemma 3.3 Let En be defined by 

1 — ^0{X) 

and let be the greatest convex minorant of 

Ho{x) + n-^/^En{x), X G [0, a]. 

Then 



^ EYn log n _ ^ flogn 



n 



[ \Mn{x) - Hn{x)} dFoix) 
J[0,a] J 

= / {Ho{x) + n-^/^En{x) - C^{x)} dFo{x) + Op 

J[0,a] ^ J 



log n 



n 



Proof. The result immediately follows from (3.20) and Lemmas 3.1 and 3.2(i). 



(3.20) 



(3.21) 



(3.22) 



□ 



We now note that the process 



X ^ . ^ r/{ , x G [0,a], 



1 - Fo{x) 



has the same distribution as the process 



X ^ Vn{x) '^^^Ho{x)+n-'/^W (j^^^) > ^ e [0,a], 



(3.23) 



where W is standard Brownian motion on M+. So, if C„ is the greatest convex minorant of the 
process 

X ^ Hoix) + n-^/^W ( ,xe[0,a], 



\l-Fo{x] 



we have: 



/ {Ho{x) + n-^^^En{x) - C^{x)} dFoix) 

J[0,a] ^ 



(3.24) 
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Theorem 3.1 Let Hq be strictly positive on [0,a], with a strictly positive continuous derivative 
on (0, a); which also has a strictly positive right limit at and a strictly positive left limit at a. 
Moreover, let Sn be defined by 

Sn{x) = Ho{x) + n-'/^W (^^^M^j - Cn{x), X e [0,a], (3.25) 

where Cn is the greatest convex minorant ofVn defined in (3.23) and let be defined by 



a 



Dn= / Snix)dFo{x). (3.26) 
Jo 

Finally, let C(0) and be defined as in Theorem 2.1. Then: 

„5/6 _ ^^^1 ^ ^(Q^ ^2^j^ n^^, 

where 

n'/'EDr, ^ E\Cm £ (^^^tII^) dHoit), (3.27) 



and 



a h!^^ ^'-^'^ 



The foUowing coroUary is immediate from the preceding. 

Corollary 3.1 Let ho be strictly positive on [0, a], with a strictly positive continuous derivative H'q 
on (0, a), which also has a strictly positive right limit at and a strictly positive left; limit at a. 
Then: 

n'/^ [Wnix) - Hn{x)] dFoix) - EDnj A iV(0,a|,J, n^oo, 

where EDn and ajj^ are defined as in Theorem 3.1. 

For the proof of Theorem 3.1 we divide the interval [0, a] into rn„ intervals /„,fc with (eqiial) 
length of order n~^/"^logn (big blocks), separated by intervals J„ ^. (/c = 2,3, . . . ,m„) with length 
of order 2n~^^^^\ogn (small blocks). The small interval Jn,\ to the left of 1^,1 has half the length 
of the other separating blocks as has the small interval Jn,m„+i to the right of In,mn- Hence, 

[0, a] = Jn,l U In,l U Jn,2 U I„,2 • • • U Jn,m„ U /n,m„ U Jn,m„+1- 

For = 2, 3, . . . , rfifi^ let Jyi ^ be the interval with the same right endpoint as Jn.k with half the 
length of Jn^k and take J„^i = Jn,i- For k = 1,2, .. . ,m„ — 1 let Jn,k+i be the interval with the 
same left endpoint as J^.fc+i with half the length of J^.fc+i aiid Jn,m„+i = Jn,m„+i- Then 

[0, a] = Jn,l U In,l U Jn,2 U Jn,2 U In,2 • • • U Jn,m„ U /n,m„ U Jn,m„+1 

where all /-intervals have the same length, of order n"~^/^logn and the J-intervals have the same 
length of (smaller) order n~^/^^/\ogn. Finally, let the interval Ln,k be defined by 

Ln,k = Jn,k U /n,fc U Jn,fc+i = Kfc, a„,fc+i), A; = 1, 2, ... , m„, yielding [0, a) = U^^^L„_fc. (3.29) 
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Note that m„ ~ an^/^/logn and see the figure below for the structure of the partition. 

Jn,l In,l Jn,2 Jn,2 In,2 Jn,3 In,mn '^n,m„+l 

1 1 1 1 1 







L. 



n,2 



The (key) locaUzation lemma below and proved in the appendix, shows that on intervals In,k 
the global convex minorant of Vn (defined in (3.23)) over [0, a] coincides with high probability with 
the restriction to I^k of the local convex minorant of the process Vn on the interval L^^k- 

Lemma 3.4 Let ho be strictly positive on [0, a], with a strictly positive continuous derivative h'g on 
(0, a), which also has a strictly positive right limit at and a strictly positive left limit at a. Then: 

(i) The probability that there exists a k, 1 < k < mn, such that the greatest convex minorant Cn 
of Vn is different on the interval Ink from the restriction to Ink of the (local) greatest convex 
minorant ofVn on Lnk, is bounded above by 

ci exp|-C2(logn)^/^| , 

for constants ci,C2 > 0, uniformly in n. 

(a) The probability that there exists a k, 1 < k < rUn, such that Cn has no change of slope in an 
interval Jnk or Jnk is also bounded by 

ciexp|-C2(logn)^/^| , 

for constants ci,C2, uniformly in n. 

For each n > 1 and 1 < A; < m„ define independent standard Brownian motions Wni , ■ ■ ■ , Wn,m„ 
and consider the processes 

X ^ Ho{x) - Hoiank) + n-'/'Wnk f " T^%^) > ^ ^ ^n.- 

\l-Fo{x} l-Fo{ank)J 

Denote the greatest convex minorants of these processes (on Lnk) by Cnk- Furthermore, define the 
processes Snk by 

Snkix) = Ho{x) - Hoiank) + n-'/^Wnk [y^^) ~ 1-FoKfc) ) " ^ ^^-^^^ 

Lemma 3.5 Assume that the conditions of Theorem 3.1 are satisfied. Moreover, letC{0) be defined 
as in Theorem 2.1 and as in Theorem 3.1. Then: 

ran 

n"'"y. I i^nk(X) - -fc'^nfe(a^)| amx) > iV(U,(7|f J, n ^ DO, 



^5/6 ^ / i^^^(^) _ ESnkix)} dFoix) A iV(0, a%J, 
k=l 



where (see (3.27)) 

1/3 



n2/3V/ ESnk{x) dFoix) -> E\Cm f dHo{t),n^oc, (3.31) 
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Proof. Let Cn = n^/^\Lnk\ ~ logn and Ink = [o-nk + n ^^"^ ^/cn.ank + n ^/^(c„ - -^/cn)]- We then 
have: 

Ink 

\ - Fo(a„fe + n-y^x) 1 - Fo(a„fc) ) "'^^ "'^ ^ J J 

• /o [a-nk + rT^/^x^ dx, 

Here we use that the (first two) deterministic terms in (3.30) drop out because of subtraction of 
the expectation. This imphes that 

n / {Snk{x) - ESnkix)} dFo{x) 

Jink 



'Ml - Fo(anfe + n-y^x) 1 - Fo(a„fe) ) '^"'^^'^^ 



i:^ J i/6w / ^o(anfc + ?i ^/^a;) i^o(anfc) \ w / nI 1 / , -1/3 \ J 



\^ 1 - Fo(anfe + n-V3a;) 1 - Fo(anfc) ^ 
where C^^ is the greatest convex minorant of the process 

X I-)- n^/^ |i7o(anifc + n~'^/^x) - Ho{ank) - n~^/^x^o(onfc)} 



yi- Fo{ank + n-^/^x) 1 - i^o(anfc) / 

Here we use that adding a hnear function to a function does not change the difference between 
this function and its greatest convex minorant. Note that the integrals on Ink only depend on 
the increments of the Brownian motion process on the corresponding disjoint intervals Lnk and 
therefore are independent. For the individual integrals we are close to the situation of Theorem 
2.3, with, for c„ oo, on [0, c„] (note that n is determined by c„, n = e'^"-) 

F^Jx) = nV3 |Fo(a„fc + n-^/^x) - Fo(a„fc)} = fo{ank)x{l + o(l)), 

Hc^{x) = r?!'^ {i/o(anfc + rT^I'^x) - i/o(anfc) - n-^/3x/io(anfc)} = ^/io(anfc)x^(l + o(l)) 
A I \ 1/3/ F(i(ank^n-^l'^x) F^iank) \ foiank)x , , 



This yields: 

n 

var ■ 



/ Snk{x) dFQ{x) j ~ (j^jt, n ^ 00, 
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uniformly in A; = 1, ... , m„, where 

(/o(anfc)/{l - Foiank)}'^)^^^ foiankf 



^2 _ 
'^nk — 



a 



777; Cr = —-, TTTTT (T 



;^o(anfc)) 



^'o(a„fc)4/3 



and cr^ is defined as in Theorem 2.1. Likewise, also with C(0) as defined in Theorem 2.1, 



n 



ESnk{x) dFQ{x) 



^ Jink 

Since the fourth moments of 



2^l^hQ{ank){hQ{ank)Hank)Y'^E\Cm 
Kiank)'/' 



/ {Snk{x) -ESnkix)} dFo{a 



are uniformly bounded by Theorem 2.3, we get for each e > 0, using Chebyshev's inequality. 



fc=i 



-1/2 



n 



/ {Snk{x) - 

J Ink 



ESnkix)} dFoix) 



> £ > ^ 0, n ^ oo. 



Using that to„^ ^ a ^/^logn and that the intervals Ink have lengths of order n ^^^logn, we 
get: 

,1/3 



m. 



' ^ KM^^ 



k=l 



k=l 



2^/3^2 r ho{tf{hoit)fo{t)y/' 



dt. 



a Jo Kitr/^ 

Since m„ = an^^^/cn, the normal convergence criterion on p. 316 of Loeve (1963) now gives: 

^5/6 ^ / i^^^(^) _ ESnkix)} dFoix) 
k=l •'^■"■k 

= m-V2 ^ !^ / - ESnkix)} dFoix) A iV (O, . 



Also note that: 



2^/^/to(anfc) {hoiank)foiank)Y'^ 
Kiank)^'^ 



''t"n f '"'n 

-n" E ^ / dFoix) ~ m-V24/2 ^ 

, /^z^ir^^n^i r '^"'hoit){hoit)foit)Y'\ /7 2/^o(t)/o(t) y/^ 



□ 



In applications of this result, used in a bootstrap approach to the computation of critical values, 
we need the following lemma, which gives a more precise expansion of the asymptotic representation 
of the expectation, given in (3.31). 
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Lemma 3.6 Assume that the conditions of Theorem 3.1 are satisfied and assume in addition that 
ho has a hounded second derivative on [0,a]. Then 



n'l^ED^ = E\Cm £ (^^^^m^) dHoit) + o (n-V^) . (3.32) 



Proof. We have: 

1/3. 



n 

c, 



- ESnk{x) dFQ{x) = 5^ ^y/^ + o\n ' j , (3.33) 



uniformly in k = 1, . . . , This is seen in the following way. 
On the intervals Ink we get: 

Fcjx) 4^^nV3 {Fo(anfc + ^"^3^) _ Foiank)} = fo{ank)x + O (n-V3(iog„)2^ ^ 

Hcnix) ^^n^/^ |i7o(a„A: + n'^^^x) - Ho{ank) - n~^/^xho{ank)^ 
= i/i;(a„fe)x2 + o(n-V3(logn)3), 

^"^^ " \ 1 - Foiank + n-V3x) 1 - Fo(a„,) / " (1 - Fo(a„,))2 + ^ 1" ^^^^"^ J " 

uniformly in = 1, . . . , to„. The relation for and Gc„ immediately follow from the mean value 
theorem, applied on the remainder term, together with the conditions of Theorem 3.1, which yield 
that /ig and /g are uniformly bounded. In the expansion of iJc„ we use the boundedness of the 

second derivative h'^. 



Combining these relations gives (3.33), and hence: 

,2/3 

— It , Cn 

fc=l " k 

1/3 



n 



/ ESnkix) dFoix) = -Y, ESnkix) dFoix) \Ink\ 

Jink ^'^ Jink 

n ^ ho{ank){2K{ank)h{ank)Y'^E\Q{S))\ l^n/^-i/3n ^3 



where, in the last line, we use again the boundedness of /q, /iq and h^, combined with the mean 
value theorem on the intervals Ink- Note that part (i) of Lemma 3.4 tells us that the probability 
that Sn is different from Snk on Ink is bounded above by 

ci exp |-C2(logn)^/^| , 

so we also have: 

n'"Y.I ESn{x)dFo{x) =E\CmY.I ( ''"^;;^;'(^^ ) '^^ffo(0 + O (n-V3(logn)3) . 
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It is clear that we get in a similar way: 

„2/3^ / ESn{x)dF^{x) 



Hence 



n'l^EDr. = n'/' g ES^{x) dFo{x) = £ [^^^^j^) dHo{t) + o (n-V^) . 



□ 



We can now prove Theorem 3.1. 
Proof of Theorem 3.1. By Lemmas 3.4 and 3.5 we have: 

„5/6 ^ / _ ESnkix)} dFoix) A iV(0, ajjj, n^oo. 

k=l •'^">' 

By part (i) of Lemma 3.4, the probability that Sn is different from S^k on 1^^ is bounded above by 

ciexp|-C2(logn 

implying that also: 



^5/6 J2 / {Snix) - ESnix)} dFoix) A iV(0, ajjj, n^oo. 

k=l ^■^1' 

For similar reasons we have: 

^5/6 ^ / _ ESnix)} dFoix) A 0, n oo, 

where we use Theorem 2.3 (this is the essence of the "big blocks, small blocks" method). The result 
now follows, since 



Dn = n^/' E / {^"(^) - ESnix)] dFoix) 

k=l •'^rik 



□ 



4 Further central limit results 

In order to derive the asymptotic distribution of the statistic [/„ defined in (1.2) and used in the 
simulations in Groeneboom and Jongbloed (2011a), we first consider the statistic 



/ {¥nix)-Fnix)} dFoix), 
J[0,a] 
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which is analogous to the statistic discussed in the preceding section, but has F„(x) — Fn{x) as 
integrand instead of IHI„(a:) — Hn{x). We have, if En again denotes the empirical process, and 
arguing as in the proof of Lemma 3.2 (i), 

F„(x) = 1 - exp {-M„(x)} = 1 - exp ^ -Ha{x) + \og\l " ' 



l-Fo(a;) J J 

= 1 - exp I -iJo(x) - I + Op (n-^) , 

uniformly for x G [0, a]. Hence, defining, as in Lemma 3.2, as the greatest convex minorant of 
the process 

we get, by Lemma 3.2, 

F„(x) - = exp{-i/„(x)} - exp |-ifo(x) - + 

= exp{-i/„(x)} |l - exp |-iJo(x) - ""^^'^'^^'^^ + ^^n(x)| | + {n^) . 

Next, replacing by i?„(Fo(x)), as in Lemma 3.3, where (-Bn) are the approximating Brownian 

bridges, we get: 

F„(x) - K{x) 

= exp {-C^ix)} { 1 - exp {-ilo(x) - ^^^^^^^^ + C^^x)^] + (i^i^) , 

where is the greatest convex minorant of the process 

, , n-^/^Bn(Fo(x)) ,^ , 
x ^ Ho{x) + — ;y " , X e [0,a]. 

Again using the results of the preceding section, it is seen that this implies that 
F„(x) - Fn{x) 

^ exp{-C„(ar)} |l - exp |-i?o(^) - n'^/'^W (3^^^) + Cn{x)^^ + Op (^-^^ , (4.34) 

where W is standard Brownian motion on [0, 00), and C„ is the greatest convex minorant of the 

process 

x^Ho{x) + n-'/'w(^-^^^^ ,xe[0,a]. 
This representation suggests to consider 

exp{-i?o(x)} I^Hoix) + n-'/^W (^^^^) - Cn(x)} 

= {1 - Foix)} !^Hoix)+n-^/^W (y^^^) - , x G [0,a]. 

We have the following result. 
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Lemma 4.1 Let /iq be strictly positive on [0,a], with a strictly positive continuous derivative Iiq 
on (0, a), which also has a strictly positive right limit at and a strictly positive left limit at a. 
Moreover, let Cn, Sn and Vn be defined as in Theorem 3.1 and let be defined by 

D^°= r Sn{x){l-Fo{x)}dFo{x), (4.35) 
Jo 



Then: 



where 



_ 2 r ( 2ho{t)foit) Y/' 



4 



and is defined as in Theorem 2.1. 

Proof. The only difference with Theorem 3.1 is that dF^it) is replaced by {1 — F^it)} dFQ{t) in 
the integral. This means that instead of EDn we get: 

and instead of we get: 

''h,{tf{h,{t)my'\, ^u^^2,, 2 r f2hoit)Mt)Y/^ 



- ' " X — — " "^'^'^^ - " X l^^J '"^^^^^ 



□ 



We have the following corollary. 

Corollary 4.1 Let Hq be strictly positive on [0, a], with a strictly positive continuous derivative h^ 
on (0, a), which also has a strictly positive right limit at and a strictly positive left; limit at a. 
Moreover, let S'^ be defined by: 

S'^{x) =¥n{x) - Fn{x), X e [0,a], 
where F^ is defined in (1.2) and let D'^ be defined by 

D'^= r S'„{x)dFo{x), (4.36) 
Jo 

Then, 

_ ED'^} A N{0,a%), n ^ oo, 
where ap^ is defined as in Lemma 4-1- 
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Proof. This is (in a sense) an application of the delta method. By (4.34) we can replace F„ — F„ 
by: 

exp{-C„(x)} |l - exp !^-Ho{x) - n^l'^W + C„(x)|| . 

We also have, using notation of the same type as in the proof of Lemma 3.5, 

Jo 

~ XI / E ^Ho{ank + n~^^^u) - Ho{ank) - Cn{ank + n~^^^u) + Cn{ank)j h{ank)du 
k=l •'^ 

y2 E {^hQiankW - Cnk{u)} fo{ank)du, 
where Cnk is the greatest convex minorant of the process 

X ^ iKiar^kW + W (i^!!'^^) , ^ e [0,c„]. 
By Brownian scaling, we get: 

E{lhQ{ankW - Cnkiu)}"^ h{ank)du 





-^o(anjk)/ 

where C is the greatest convex minorant of x W{x) + x^, a; G M. So we find: 

£e{Ho{x) - Cn{x)f dFoix) ~ n-^/'EC{Of£ {lh',{t)f' { l-F^it) ^ ^^'^^^ 
We also have: 

£ E I^Hoix) + n~y^W (^l^^j^) - ^n(x)|' dFoix) 

~ n-'/'EC{or£ mit)f' (r^^)'^' ^'-^^^ 

Hence, by (4.37) and (4.38), 
J\xp{-Cnix)} |l - exp |-go(:^) - n-'/^W +C„(x)|| dFo{t) 

= [{^- |l - exp |-iJo(x) - n-^/''W (j^^) + } ^^^0(0 + Op 

where we also use the Cauchy-Schwarz inequality in the first equality. 
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For the expectation we get similarly 

'logn 



£ E [Wnix] - Fn{x)] dFoix) = ED^° + O 



n 



where is defined by (4.35). This is seen in the following way. Since we assume that Fq^o) < 1, 
we have by by Chernoff 's theorem (as in the proof of Lemma 3.2), 

P (1 - ¥n{a) < i{l - Fo(a)}} < 6""^ 

for a c > 0, and hence, defining the event An by 

An = {l-¥n{a)>^{l-Foia)}}, 

we get 

£ E\¥n{x) - Fn{x)j dFo{x) = £ E\¥n{x) - Fn{x)^ 1a„ dFoix) + O (e-"^) 
= £ E |e--^"(^) - e-^"(^)} 1a„ dFo{x) + O {e''"') 
= £e{1- Fnix)} |e-{^"W-H"W} - l} dFoix) + O (e""^) 
= £{^- M^)} E |e-{^"(^)-="(^)} - l} 1a„ dFoix) + O [n'^) 
= £{^- Foi^)} E iyHoix) + n-^/'^W ( J^^) - ^n(x) | dFoix) + O 



logn 



£;£>f° + O 



log n 



n 

The result now follows from Lemma 4.1. □ 

Similarly as in Lemma 3.6, we have the following expansion for the expectation in Corollary 
4.1. The proof proceeds along the same lines as the proof of 3.6 and is therefore omitted. 

Lemma 4.2 Assume that the conditions of Corollary 4-1 o-f^ satisfied and assume in addition that 
ho has a hounded second derivative on [0,a]. Then 

r^-l-ED'^ = E\CiO)\ £ (^^^^111^) dFoit) + a (n-V^) . (4.39) 

The preceding results finally yield the following theorems. 

Theorem 4.1 Let be defined as in Theorem 3.1 and let the conditions of Theorem 3.1 he 
satisfied. Then: 

n^/e {m„(x-) - Hnix)] (Mnix) - EDn^ A Ar(0,(7|,J, n ^ oo, 
where ED^ and a^^ are defined as in Theorem 3.1. 
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Theorem 4.2 Let the conditions of Lemma 4-1 be satisfied and let cr|.p and C(0) he defined as in 
Lemma 4-1 ■ Then: 



„5/6 



and 

We only prove Theorem 4.1, since the proof of Theorem 4.2 proceeds along similar lines. 
Proof of Theorem 4.1. Using Lemma 3.1 again, we get: 



j\nn{x-)-Hn{x)] d¥n{x) ^ j\Ho{x)+n-^'^Wn {j^^^yj -Cn{x)] dFn{x)+Op 



logn 
n 



where C„ is the greatest convex minorant of Vn which is defined as in (3.23) with Wn replacing 
W. The process W„ is distributed as standard Brownian motion on [0,a] and Wn o {Fq/{1 — Fq)) 
is given by 

/ Fo(.) \_ Bn{F,{x)) 

where Bn is coupled to the empirical process as in Lemma 3.3. 
We only have to show 

{Vn{x) - Cn{x)] d (F„ - Fo) {x) = Op (n^^/ej ^ (4 40^ 



since we then have: 



£ {Hn(x-) - Hn{x)] d¥n{x) = £ {Vn{x) - Cn{x)} dF„(x) + Op 



logn 



n 



logn 



r {Vnix) - Cnix)} dFo{x) + T {K(x) - C„(x)} d (F„ - Fq) (x) + Op 

^0 
= £ {Vnix) - Cnix)} dFoix) + Op (n-^/^) . 

To show that this relation holds, we follow a method which is somewhat similar to the method 
used in Kulikov and Lopuhaa (2008) (but uses the Brownian motion representation instead of the 
empirical process and does not bring the derivative of the greatest convex minorant into play). 
The ^?-variation of a function / on the interval I = [0, o] is defined by 

I) = sup 1^ Ifi^i) - fixi-i)f : xo = < xi < ■ ■ ■ < Xm = . 
The variation norm of / on I is defined by 

ll/ll[p] = ^p(/;^)'/^ + sup|/(x)|. 
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We have, by Theorem II.3.27 in Dudley and Norvaisa (1999), for p,q> and 1/p + l/g > 1: 



/ 



{Vn{x) - Cn{x)} d{¥n-Fo){x) 



< C \\Vn — Cn||[p]||F„ — Fo\\[q], 



(4.41) 



for a constant c > 0. Moreover, by Theorems 1.6.1 and 1.6.2 in Dudley and Norvaisa (1999), and 
Theorem 3.2 in Qian (1998), we have: 



ll^n-i^0||[5] 



' Op (n^i-")/'?) ,qe[l,2), 

0^ [n-y^^/L{Ln))) , g = 2, 
. Op (n-V2) , ^ > 2, 



(4.42) 



where Ln = 1 V log n. 

Let Ti, . . . , Tm be the points of jump of the derivative c„ of Cn on [0, a], and let tq = 0, Tm+i = o- 
The function C„ is linear on the intervals [ri,ri_|_i], and Vn behaves on such an interval as an 
excursion above its greatest convex minorant C„, with the same values as Vn at the endpoints of 
the interval. Hence we have, for p > 2, by Lemma 4 of Huang and Dudley (2001), 



m+l 



m+1 



i^piVn - C„; [0, a]) < 2P-1 J] Up{Vn - C„; r,]) = 2^"^ z/p(F„; r,]). 



i=l 



where, using the fact that the linear part drops out in taking the comparison with the greatest 
convex minorant. 



Foix) 1 - Fo{Ti_i) J ' Ti - Ti-i \l-Fo{n) l-Fo(r,_i) 

X — T' 1 

+ Ho{x) - Ho{Ti-i) — {Ho{Ti) - Ho{Ti-i)} , X e [Ti_i, Tj" 

n — Ti-i 

By part (ii) of Lemma 3.4 we have: 

^max(ri — Tj-i) = O {^"^^^ logn^ . 
Let Ui be the midpoint of the interval [Ti_i,ri] and let Jho by defined by 

fHo{x) = Hq{x) - Ho{Ti-i) — {Ho{Ti) - Ho{Ti-i)} , X e [Tj_i, Tj] . 



n - Ti-i 



Then 



ho{Ui){x - Ti-l}{Ti - X}{1 + Op(l)}, 



where x ^-^ {x — Ti_i}{rj — x} is increasing on [Ti-i,Ui] and decreasing on [ui,Ti], and 

if Ho-, [ri-i,n])) ~ 2^-Ph'o{uir - T,_i)r {n - ui)r , 
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(see, e.g., (3.4) of Huang and Dudley (2001)). Hence, for any p > 2, 

m+l 

J2 ^p(M; [n-i,Ti]) 

i=l 

m+l m+l 

~ 2'-p h'oiuif {ui - n-iV {n - = 2'-'p ^ol^iF {^i - n-if" 

1=1 i=l 

m+l 

i=l 

~ mpc{ri - T,_i}2f-i h'^iuf du = (n-(2p-i)/3(logn)(2f-i)/2) . 

Note that the Op-term becomes Op (n~^(logn)'^/^) for p = 2. 
For the Brownian part 



l-Fo(a;) l-Fo(ri_i)y - r^-i \1 - Foin) 1 - Foin-i) 



we find, for ^? > 2, 

m+l 

-p/2 



XI [Tj-l,rj]) = Op 

1=1 



by the fact that almost all Brownian motion paths are Holder continuous of any order < 1/2. 
So we find: 

UK - C„||[p] = Op (n-V2(iogn)(2^'-i)/(2^)) , (4.43) 
for any p > 2. Thus (4.41), (4.42) and (4.43) imply 

/ {Vnix) - Cnix)} d (F„ - Fo) (x) = Op (n-i+^) , 

J[0,a] 

for arbitrarily small £ > 0. □ 

We end this section with a result for the situation that the hazard is nondecreasing, but not 
strictly nondecreasing. 

Theorem 4.3 Let Fn and F„ he defined as in Theorem 4-^ and let (again) 

Un= [ {¥n{x-)-Fn{x)}dFn{x). 
J[0,a] 

Let U he given by 

U = £ll- F„(.)} [w (^^) - C(x)} 
where W is standard Brownian motion on [0, oo) and C is the greatest convex minorant of 

x^w(-^^^^^],xe[0,a]. (4.44) 



l-Fo{x) 

Suppose that the underlying hazard Hq is constant on [0, a] . Then: 

n^^^Un U, oo. 
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Proof. The proof follows lines that are familiar by now. We first consider 

Un= [ {¥nix-)-Fn{x)}dFoix). 
J[0,a] 

By (4.34) we can replace F„ — F„ by: 

exp{-C„(x)} |l - exp !^-Ho{x) - n-'/^W (j^^^) - Cn(x)|| , 

where Cn is the greatest convex minorant of the process 

X ^ Hoix) + n-'/'W ( J^^) > ^ e [0, a], 

with a remainder term of order Op ( (log n)/n). Using the delta method as in the proof of Corollary 
4.1, we can replace this (apart from a remainder term of order Op{n~^)) by: 

n-V2|i _ Foix)} (l^^) - C(x)} , X G [0,a]. 

where C is the greatest convex minorant of the process (4.44), and where we use that Hq is linear 
on [0, a]. The statement for C/„ now follows by an application of Dudley and Norvaisa (1999), as 
in the proof of Theorem 4.1. □ 



Remark 4.1 Note that the rate of convergence drops from n^/^ to n-^/^ in Theorem 4.3, and that 
the limiting distribution is not normal. We get a limit behavior that can be analyzed using the 
methods of Groeneboom (1983), where the concave majorant of Brownian motion without drift is 
characterized via a Poisson process of jump locations and Brownian excursions. 



Appendix 

Proof of Lemma 2.1. Let u > 0. Then, for x > u: 

V{x) = W{x) + {x-uf +2u(x-u)+u'^ > W{x) + (x-uf+u^ = W{u)+u^+W(x)-W{u) + {x-uf. 
Hence, 

P I miny(x) < I < P ( min W{u) + u'^ + W{x) - W{u) + {x-u)^ <0] = 



P W{u) + u^ + mmW{x) - W{u) + {x-uy <0] < 

\ x>u J 

P {W{u) < -^u^) +¥(uimW{x) - W{u) + {x - uf < -^u^^ 
The process 

X W{x) - W{u) + {x- uf, x>u 

behaves in the same way as the process 1 1— t- V{t), t > 0, but starts in x instead of 0. By Corollary 
2.1 in Groeneboom and Temme (2010) we have that for all z > 0, 

P jmin V{t) < -z| ~ 2 • 3"^/^ exp |-8z^/V^} , ^ ^ oo, (4.45) 



26 



implying that there exist positive constants ci and C2 such that for all u > 

P (mmW{x) - W{u) + {x - uf < - W] < ciexp{-C2U^} . 

\x>u J 

We also have for all u > 

P {W{u) < -lu'} = P {Win)/V^ < -i^^/n < 

implying that there exist positive constants C3 and C4 such that for all tt > 

F{W{u) < -^u^} < C3exp{-C4U^} . 

Combining these upper bounds with the fact that the process V running to the left from zero 
behaves in the same way as the process V running to the right from zero, the statement of the 
lemma follows. □ 



Proof of Lemma 3.4. (i). The interval In,k is bounded on the left by the interval J„ ^ and on the 

right by the interval Jn,fe+i- The intervals Jn,k and Jn,k+i both have length of order n~^/^-^logn. 
If the greatest convex minorant C„ of Vn on [0, a] has changes of slope in the intervals J^^k and 
Jn,k+i, the greatest convex minorant of Vn on [0,a], restricted to the interval In,k, coincides with 
the greatest convex minorant C^k of Vn on -L^^fc, restricted to the interval In,k- So we have to find 
bounds for the probability that the greatest convex minorant of Vn on [0, a] has no changes of slope 
in Jnk or Jn,k+i- To do this, we follow the method used in Groeneboom and Wellner (1992), p. 
96. 

Let ttnk and bnk be the left and right endpoints of Jn,k+i, respectively, and let Unk be its 
midpoint. If 

Cn{ank) < ho{Unk) < Cn{bnk), (4.46) 

where c„ is the left-continuous slope of (7„, then Cn has a change of slope in the interval Jn,k+i- 
Note that for x > bnk, using the assumed smoothness of Hq, and inf[o,a] ^oi^) = 2k > 0, 

Vnix)-Vn{Unk) > n'^^ (j^^^) " ^ ( ^ ^^^(ti) ) } + ^o(^nfc) - Unk) + <X - Unk?- 

(4.47) 

Now consider the event that 

Cn{bnk) < hoiUnk), (4.48) 

and let Tnk be the first point of jump of c„ to the right of bnk- Then 

Cn{x) < hoiUnk), X < Tnk, 

and hence 

/•Tnk 

Vn{Tnk) - Vn{x) < Cn{Tnk) - Cn{x) = I Cn{y) dy < ho{Unk){rnk - x), X < Tnk- 

J X 

Using (4.47) and stationarity of Brownian Motion, this means that the probability of (4.48) is 
bounded above by 

^{yn{Tnk) - VniUnk) < ho{Unk){Tnk - Unk)} < P {3a; > bnk ■ Vn{x) - Vn{Unk) < ho{Unk){x - Unk)} 

. r : {w{^yw (^) } . - 
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We will see that this probability will become exponentially small. To this end, define the following 
covering of [6„fe,a] 

Knkj =^ [tnkj,tnk,j+i\ =^ + jn~^^^, h^k + (j + l)rr^/^ = [bnk, a] 

for < j < [n^/^(a — bnk)\ (where the right end point of the last interval is taken to be a). Then 
the probability in (4.49) can be bounded above by 

E ' {^^ ^ : -"^ (i^ - 1^) } <- - ■ 

Denoting the probabilities in this sum by Pnkj , we get 

Pnkj < P < sup W ( - r^^iTTT^ ) ^ I^V^itnkj - Unkf \ 

<¥< sup W{z) > nVnitnkj - u^k? \ ■ 

[o<z<Fo(t„fc,,+i)/(l-Fo(t„fc,,+i))-Fo(w„fc)/(l-FoKfc)) J 

Since tnk,j+i & [^nfcifl] for all j under consideration, 

^ Fo{tnk,j+l) Fo{Unk) {Fo(tnk,j+l) - Fo{Unk)) _ x 

- 1 - Fo(t„,,,+i) 1-FoM- {l^Foia))^ S ^[t^k,j+i ^nfcj 

for some < A < oo, we obtain, for a standard normal random variable Z 



Pnkj < IP < sup W{z) > nVnitnkj - Unk) } = ^<\Z\> 



u r f I I \J ^ V TIKI Tlhj J i — I I I — J I 

0<z<\{tr,k,i+i-Uuk) J I, \J Ktnk,i+1 - Unk) ] 

< P||Z| > KA/n(t„fej - -U^jfc)^/^! < ]^ClC£>{-\nR^{tnkj-Unkf}- 

Using that tnkj - Unk = Kk - Unk + jn~'^/^, and bnk - Unk ~ ln~^/^^/logn, we get 

Lni/:^(a-6„fe)J 

Pnkj < exp (|(logn)3/2 + I ^ ^ < pexp |-p'(logn)3/2| 

for some p, p' > 0. Combining this with (4.49) and (4.50), this bounds the probability of (4.48) 
from above. Since a similar bound holds for the probability of the event c„(a„jk) > ho{unk), the 
probability that (4.46) does not hold for a specific k, is bounded by a bound of the same structure. 
Moreover, since this upper bound does not depend on k and nin ~ an^/^/ log n, the probability 
that there exists a 1 < k < m„ for which (4.46) does not hold satisfies the same bound (with slight 
change in p and p'), this proves (i). Part (ii) is an immediate consequence of (i). □ 
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